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SUMMARY 


An  account  is  given  herein  of  an  investigation  conducted  to  deter¬ 
mine  the  feasibility  of  artificially  thickening  a  turbulent  boundary 
layer  on  a  flat  plate.  A  description  is  given  of  several  methods  used 
to  thicken  artificially  the  boundary  layer.  It  is  shown  that  it  is  pos¬ 
sible  to  do  substantial  thickening  and  obtain  a  fully  developed  turbu¬ 
lent  boundary  layer,  which  is  free  from  any  distortions  introduced  by  the 
thickening  process,  and,  as  such,  is  a  suitable  medium  for  fundamental 
research.  Measurements  of  mean  velocity,  spectrum  of  u-f luctuation,  and 
intensity  of  u-fluctuation  which  served  as  criterions  for  determining 
the  state  of  the  layer  are  presented.  Some  features  of  the  fully  devel¬ 
oped  turbulent  boundary  layer  with  zero  pressure  gradient  are  discussed. 
The  mean  flow  in  the  boundary  layer  is  compared  with  the  logarithmic 
laws  of  mean- velocity  distribution  derived  for  pipes  and  channels,  and 
the  existing  theories  dealing  with  the  shape  of  the  spectrum  in  isotropic 
turbulence  are  applied  to  the  spectrum,  measured  in  the  outer  portion  of 
the  boundary  layer.  It  is  also  shown  that  the  turbulence  in  the  outer 
portion  of  the  boundary  layer  has  an  intermittent  character  similar  to 
that  found  for  wakes  and  jets. 


INTRODUCTION 


The  lack  of  a  satisfactory  theory  for  turbulent  shear  flow  and  the 
development  of  the  hot-wire  anemometer  as  a  powerful  tool  for  turbulence 
research  have  directed  attention  to  the  quantitative  measurement  of  the 
characteristics  of  turbulence  in  shear  flow  as  a  means  of  obtaining 
information  on  which  to  base  a  sound  theory.  The  recent  concept  of 
local  isotropy  proposed  by  Kolmogoroff  (references  1  and  2);  that  is, 
that,  whatever  the  nature  of  the  large-scale  motions,  the  small-scale 
turbulence  tends  to  be  isotropic,  has  intensified  the  experimental 
program  with  the  hope  that  the  existing  theories  of  isotropic  turbulence 
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may  be  useful.  While  a  considerable  amount  of  work  has  been  done  in 
various  flow  fields  (references  4,  5 ,  and  6),  attention  is  here 
^entered  on  the  boundary  layer,  not  only  because  of  its  practical  impor-  jr 

stance,  but  to  provide  information  which  can  be  compared  with  that  obtained 
in  other  types  of  shear  flow  and  thus  show  differences  and  likenesses 
and  in  particular  show  whether  certain  phenomena  are  universal.  It  is  V 

felt  that  the  results  of  the  measurements  can  best  be  compared  with  the 
theory  of  isotropic  turbulence  if  such  complicating  factors  as  pressure 
gradient,  curvature,  surface  roughness,  and  compressibility  are  elimi¬ 
nated.  Furthermore,  the  turbulence  should  be  fully  developed  and  the 
Reynolds  number  should  be  sufficiently  high  to  bring  the  domain  of  local 
isotropy  within  the  range  of  measurement.  In  addition  the  boundary  layer 
should  be  as  thick  as  possible  in  order  to  facilitate  the  use  of  multi¬ 
wire  hot-wire  probes  and .minimize  errors  due  to  wire  length. 

These  requirements  mean  that  the  surface  should  be  smooth  and  flat 
and  that  high  Reynolds  number  should  be  obtained  by  thickness  rather 
than  by  high  speed.  Such  requirements  are  difficult  to  meet,'  except  in 
wind  tunnels  especially  constructed  for  the  purpose.  The  use  of  one 
wall  of  a  conventional  wind  tunnel  does  not  appear  to  be  a  satisfactory 
solution  because  of  the  likelihood  of  extraneous  effects.  The  problem 
of  obtaining  a  thick  boundary  layer  was  encountered  at  the  National 
Bureau  of  Standards  when  a  turbulent-boundary- layer  investigation  using 
a  flat  plate  was  undertaken.  The  test  section  of  the  tunnel  was  19  feet 
long,  but  it  was  not  feasible  to  use  a  plate  longer  than  12  feet.  It 
was  felt  that  the  need  for  a  thicker  boundary  layer  than  one  which  would  to 

normally  develop  on  12  feet  of  surface  was  sufficiently  urgent  here,  and 
possibly  in  other  laboratories  with  even  less  length  available,  to  warrant 
an  investigation  of  the  possibilities  of  artificially  thickening  the 
layer.  This  involved  not  only  methods  pf  thickening,  but  also  setting 
up  criterions  by  which  to  establish  the  identity  of  a  fully  developed  . 
layer  over  a  smooth  surface. 

In  determining  the  feasibility  and  validity  of  artificial  thickening 
a  number  of  measurements  had  to  be  made,  including  mean  velocities,  inten¬ 
sity  of  u-fluctuations,  and  the  spectrum  of  u-fluctuations.  As  a  con¬ 
sequence  it  is  possible  to  show  how  such  quantities  are  affected  by 
thickening  devices  and  what  values  they  assume  when  the  layer  is  no 
longer  influenced  by  its  beginnings.  It  is  obviously  very  necessary 
to  know  whether  there  exists  a  "flat-plate  boundary  layer"  of  universal 
character  depending  only  on  Reynolds  number.  This  is  shown  to  be  the 
case  insofar  as  the  observed  quantities  provide  a  proper  sampling. 

The  investigation  was  conducted  under  the  sponsorship  and  with 
the  financial  assistance  of  the  National  Advisory  Committee  for 

Aeronautics.  * 
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SYMBOLS 

x  distance  along  surface  from  leading  edge  of  flat  plate 

y  distance  normal  to  surface,  measured  from  surface  of 

flat  plate 

U  mean  velocity  in  boundary  layer 

mean  velocity  in  free  stream,  just  outside  boundary  layer 

u,  v,  w  x,  y,  and  z  components  of  instantaneous  turbulent- velocity 

fluctuations 

“2  ~~2 

u,  v,  w  mean-square  values  of  u,  v,  and  w 

u*  root-mean- square  value  of  u 

p  density  of  air 

V  kinematic  viscosity  of  air 

P  static  pressure 

PQ  static  pressure  at  x  =  10. 5  feet 

q^  free- stream  dynamic  pressure 

q0  free- stream  dynamic  pressure  at  x  =  10.5  feet 

qs  dynamic  pressure  in  boundary  layer  as  measured  with  a 

surface  tube 

T  shearing  stress 

t0  shearing  stress  at  surface 
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■boundary-layer  thickness 

boundary- layer  displacement  thickness 


boundary- layer  momentum  thickness 


Reynolds  number  based  on  boundary- layer  thickness 

Reynolds  number  based  on  boundary- layer  momentum 
thickness 


longitudinal  space  correlation  coefficient  of  u-fluctuation 


longitudinal  scale  of  u-fluctuation 
microscale  in  isotropic  turbulence 

total  dissipation  of  turbulent  energy  per  unit  mass 
per  second 

time 

frequency,  cycles  per  second 
one- dimensional  wave  number  m 
three-dimensional  wave  number 

normalized  one-dimensional  spectral  function  associated 
with  frequency  n  .  • 
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Ex(kx)  one-dimensional  spectral  function  associated  with  kx 

E(k)  three-dimensional  spectral  function  associated  with  k 

APPARATUS  AND  PROCEDURE 
General  Arrangement 


'  The  investigation  was  conducted  in  the  4^  -  foot  wind  tunnel  at  the 

National  Bureau  of  Standards,  using  a  flat  plate  12  feet  long,  mounted 
vertically  along  the  center  line  of  the  test  section,  the  leading  edge 
being  6  feet  from  the  upstream  end  of  the  test  section.  The  plate  was 
l/4-inch  aluminum  with  a  smooth  surface  and  a  sharp  leading  edge  formed 
by  a  smooth  symmetrical  taper  beginning  4  inches  from  the  leading  edge. 
The  plate  was  fastened  along  the  horizontal  edges  to  one  side  of  3- inch 
channel  irons  anchored  to  the  upper  and  lower  tunnel  walls,  and  so 
divided  the  test  section  into  nearly  equal  passages.  It  was  considered 
desirable  to  conduct  the  investigation  in  a  tunnel  of  low  turbulence  so 
that  the  turbulence  of  the  free  stream  would  be  negligible  compared  with 
that  in  the  boundary  layers.  The  turbulence  level  of  the  tunnel  was 
0.02  percent  at  30  feet  per  second  and  0.04  percent  at  100  feet  per 
second.  An  elevation  view  of  the  tunnel  is  shown  in  figure  1.  More 
details  as  to  the  arrangement  of  damping  screens  used  in  obtaining  the 
low  turbulence  are  given  in  reference  7* 

In  order  to  obtain  a  condition  of  zero  pressure  gradient,  a  flexible 
second  wall  of  sheet  aluminum  was  mounted  on  the  tunnel  wall  facing  the 
working  side  of  the  plate.  It  was  attached  to  the  ends  of  screws 
threaded  into  the  tunnel  wall,  and  thus  provided  an  adjustable  section 
to  make  the  passage  between  the  plate  and  the  wall  sufficiently  diverging 
to  offset  the  pressure  drop  resulting  from  boundary- layer  growth. 

Another  important  consideration  was  proper  angle  of  attack  to 
afford  smooth  flow  entry  at  the  leading  edge  on  the  working  side  of  the 
plate.  This  was  accomplished  by  displacing  the  upstream  edge  of  the 
plate  and  by  partial  blocking  of  the  passage  on  the  working  side.  Since 
the  obstructions  causing  the  blocking  had  to  be  placed  farther  downstream 
than  any  usable  portion  of  the  plate,  the  passage  was  extended  4  feet 
by  joining  the  trailing  edge  of  the  plate  with  a  section  of  plywood. 

Under  these  conditions  the  pressure  distribution  along  the  plate 
was  found  to  be  as  shown  in  figure  2(a).  From  2  feet  on  downstream 
the  pressure  is  seen  to  be  uniform  on  the  average  with  Variations  of 
about  1  percent  of  q0  about  the  mean.  These  variations  could  not  be 
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removed  by  the  adjustable  wall  because  they  were  associated  with  a 
slight  waviness  of  the  surface  of  the  plate. 

Figure  2(b)  shows  the  transition  region  at  a  speed  of  108  feet 
per  second  for  smooth  entry  conditions  at  the  leading  edge.  The  begin¬ 
ning  of  transition  corresponds  to  a  Reynolds  number  (based  on  length  of 
laminar  flow)  of  3,175*000,  and  here  transition  results  from  the  growth 
of  more  or  less  regular  laminar-boundary- layer  oscillations. 


Instrumentation  and  Techniques 

Traversing  mechanism.-  In  order  to  traverse  in  a  horizontal  plane 
normal  to  the  surface  the  unit  shown  in  figure  3  was  used.  A  micrometer 
screw  with  0. 5-millimeter  pitch  is  driven  through  a  set  of  beveled  gears 
and  a  shaft  by  a  selsyn  motor  remotely  controlled.  To  provide  for  ease 
of  alinement,  the  pressure  probes  and  hot-wire  holders  are  attached  so 
that  they  can  move  fore  and  aft,  and  rotate  in  a  vertical  plane.  By 
means  of  the  bracket  shown  in  the  figure  the  unit  is  bolted  to  a  stream¬ 
lined  strut  so  that  the  end  opposite  the  gears  is  forced  against  the 
surface.  Channel  irons  anchored  to  the  floor  and  ceiling  of  the  tunnel 
served  as  guide  rails  for  the  strut,  and  traversing  longitudinally  was 
done  by  moving  the  strut  fore  and  aft  by  means  of  a  chain  and  sprocket 
manually  operated  from  outside  the  tunnel. 

Measurement  of  pressure  distribution,  transition,  and  mean 
velocity.-  Measurements  of  pressure  distribution  and  transition  were 
made  with  the  pitot-static  head  shown  in  figure  4.  The  impact  tube 
makes  contact  with  the  surface,  and  the  static  tube  is  opposite  the 
point  of  contact  l/4  inch  away  from  the  surface.  The  static  and  impact 
tube  are  made  of  0.04- inch  nickel  tubing  with  0.003-inch  wall  thickness. 
The  static  tube  is  made  according  to  the  conventional  design  for  such 
a  tube,  and  the  impact  tube  is  flattened  at  the  end  to  form  a  rectan¬ 
gular  opening  0.014  inch  wide.  The  static  tube  was  used  for  the  meas¬ 
urement  of  pressure  distribution  and  both  impact  and  static  tubes  for 
the  measurement  of  transition.  The  head  was  attached  by  a  long  arm 
(containing  the  pressure  leads)  to  the  strut  and  could  be  moved  back 
and  forth  in  the  x-direction  while  the  forklike  guides  maintained  proper 
relation  to  the ‘surface.  The  determination  of  the  transition  point  with 
the  surface  tube  depends  on  the  increase  in  impact  pressure  close  to 
the  surface  when  the  layer  becomes  turbulent.  The  abruptness  of  the 
increase  is  shown  in  figure  2(b).  The  beginning  and  ending  of  the 
transition  region  coincide  with  the  beginning  and  ending  of  the  rise. 

This  was  checked  with  the  hot-wire  head  shown  in  figure  4.  This  sled¬ 
like  head  was  attached  to  the  strut  so  that  it  could  be  moved  fore  and 
aft  like  the  pitot-static  head.  The  hot-wire  was  at  a  fixed  distance 
of  0.03  inch  from  the  surface,  and  visual  observation  of  the  fluctuations 
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on  an  oscilloscope  for  different  positions  along  the  surface  confirmed 
the  indications  of  the  surface  tube.  Measurements  of  mean  velocity 
were  made  with  a  pitot- static  combination  similar  to  the  one  described 
above  but  without  the  guides.  It  is  shown  with  the  traversing  unit  in 
figure  3.  No  correction  was  made  for  the  effect  of  turbulence  on  the 
readings . 


Measurement  of  intensity  and  spectrum  of  u-fluctuation.-  A  descrip¬ 
tion  of  the  method  and  the  hot-wire  equipment  used  for  measuring  the 
intensity  of  u- fluctuation  is  given  in  detail  in  reference  8.  Platinum 
wire  0.0001  inch  in  diameter  and  0.5  millimeter  long  was  used.  This  was 
wire  drawn  by  the  Wollaston  process ,  and  the  method  of  making  a  hot-wire 
probe  was  to  etch  away  the  silver  coating  and  then  soft- solder  the  wire 
to  the  tips  of  fine  sewing  needles.  The  time  constant  of  the  wire  was 
approximately  0.25  millisecond  and  was  determined  by  the  square-wave 
method  described  by  Kovasznay  (reference  9)*  Compensation  for  the  time 
lag  was  made  by  a  resistance-capacitance  network  in  the  amplifier.  The 
response  of  the  uncompensated  amplifier  was  flat  from  10  to  5000  cycles 
per  second  and  fell  off  to  about  50  percent  at  10,000  cycles  per  second. 
The  mean- square  output  was  read  on  a  microammeter  working  out  of  a 
thermocouple  circuit. 


In  order  to  measure  the  spectrum  of  u-fluctuation  the  output  cir¬ 
cuit  of  the  hot-wire  amplifier  was  bypassed  and  the  signal  fed  into  an 
Electrical  Research  Products,  Inc.,  sound  analyzer  with  a  frequency 
range  from  10  to  10,000  cycles  per  second  and  a  choice  of  fixed  band 
widths  of  nominally  5  or  50  cycles  per  second.  The  rectifier  type  of 
output  circuit  of  the  analyzer  was  replaced  by  an  additional  stage  of 
amplification  and  a  thermocouple  and  meter.  The  purpose  of  this  was 
to  read  mean- square  currents  and  also  to  obtain  a  highly  damped  indi¬ 
cator.  In  reading  the  output,  averages  were  taken  over  a  period  of 
1  minute.  The  mean- square  voltage  corresponding  to  the  compensated 
hot-wire  signal  passed  by  the  band  was  obtained  by  first  reading  the 
hot-wire  signal  on  the  output  meter  and  then  matching  this  signal  with 
a  known  sinusoidal  signal  fed  into  the  system  from  a  General  Radio 
micro volter.  The  known  signal  had  the  frequency  corresponding  to  the 
center  of  the  band.  The  mean- square  hot-wire  voltage  per  unit  frequency 
was  then  obtained  by  dividing  by  the  effective  band  width.  The  effective 
band  width  was  defined  as  the  rectangular  band  width  having  the  same 
area  as  the  experimentally  determined  band  shape.  These  were 


Effective 

(cps) 


Nominal 

(cps) 


7-07 

57-9 


5 

50 
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It  was  found  to  be  advantageous  to  use  the  7 • 07  band  from  10  to 
1000  cycles  per  second  and  the  57-9  band  from  1000  to  10,000  cycles 
per  second. 


Artificial  Thickening  of  Boundary  Layer 

It  is  at  once  obvious  that  objects  protruding  outward  from  a 
surface  create  more  drag  than  the  surface  that  they  cover,  and  hence 
offer  a  means  of  producing  a  layer  of  de-energized  flow  in  a  compara¬ 
tively  short  distance  in  the  streamwise  direction.  The  energy  of  the 
main  flow  is  diverted  into  eddy  motion  in  wakes,  and,  while  the  flow  is 
turbulent,  it  does  not  have  the  same  dynamic  structure  as  a  fully  devel¬ 
oped  turbulent  boundary  layer  over  a  smooth  surface.  The  question  is, 
after  leaving  drag-producing  objects,  what'  length  of  smooth  surface  is 
needed  for  the  layer  to  become  the  same  as  a  smooth- surface  generated 
layer?  It  would  seem  that  the  length  would  be  less  than  that  required 
to  generate  a  layer  of  the  same  thickness,  but  this  is  by  no  means  a 
foregone  conclusion. 

Since  the  ultimate  aim  of  the  work  was  to  produce  a  layer  approx¬ 
imately  3  inches  thick,  this  thickness  having  been  chosen  as  sufficient 
for  the  convenient  use  of  hot-wire  probes,  no  effort  was  made  to  find 
the  maximum  thickness  that  could  be  produced  in  a  given  length.  Emphasis 
was  placed  rather  on  simplicity  of  method  and  adequacy  of  tests. 

The  turbulent  boundary  layer  along  the  plate  was  first  obtained  by 
the  conventional  method  of  free  transition,  that  is,  from  the  instability 
of  laminar-boundary-layer  oscillations.  Then,  various  devices  such  as 
rods,  screens,  and  distributed  roughness  were  used  on  the  laminar  layer 
to  force  transition  and  de-energize  the  flow.  Attention  was  given  to 
maintaining  the  condition  of  zero  pressure  gradient  for  the  modified 
layers  by  adjusting  the  "false"  wall.  The  first  attempt  at  thickening 
the  layer  was  to  use  rods,  which  were  0.04,  0.11,  0.l8,  and  0.25  inch 
in  diameter.  These  were  placed  singly  in  the  laminar  layer  and  in  con¬ 
tact  with  the  surface  at  4  feet  from  the  leading  edge.  As  will  be  seen, 
these  were  not  satisfactory  because  of  the  length  required  for  their 
distorting  effect  to  disappear.  The  second  method  consisted  of  placing 
wire  screens  of  various  mesh  sizes  at  different  angles  to  the  plate  and 
close  to  the  leading  edge,  so  that  the  screen  and  plate  formed  a  narrow 
wedge  opening  into  the  wind.  A  24-mesh  screen  of  0.0075- inch  wire 
diameter,  which  was  1  foot  wide,  making  an  angle  of  0.5°  with  the  plate, 
and  extending  from  floor  to  ceiling,  was  a  satisfactory  thickening 
device,  but  it  was  not  considered  practical  because  of  the  possibility 
of  accidental  shifting.  The  third  and  final  method  consisted  of 
cementing  sandpaper  to  the  surface  beginning  at  the  leading  edge  and 
extending  for  distances  up  to  2  feet  downstream.  Various  amounts  of 
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sand  roughness  were  tried  -until  the  layer  was  of  the  desired  thickness. 
The  final  arrangement  of  the  roughness  was  No.  l6  floor  sanding  paper, 
cemented  to  the  plate  from  the  leading  edge  to  2  feet  downstream,  and 
extending  from  floor  to  ceiling. 

Initially  the  one -dimensional  spectrum  of  u- fluctuations  was  thought 
to  "be  the  best  over-all  test  to  be  applied  to  the  resulting  layers. 
Accordingly,  for  a  basis  of  comparison,  considerable  attention  was  given 
to  the  spectrum  of  the  fully  developed  layer  resulting  from  free  tran¬ 
sition.  It  was  found  later  that  the  mean- velocity  profile  was  a  still 
more  sensitive  test  of  the  condition  of  the  layer.  Finally,  the  dis¬ 
tribution  of  u'  across  the  layer  was  added  to  the  list  of  criterions. 
These  three,  compared  in  different  parts  of  the  layers  and  compared  with 
other  available  data,  served  as  the  sample  phenomena  to  determine  the 
state  of  the  layer. 


RESULTS  AND  DISCUSSION 
Mean -Velocity  Measurements 

With  free  transition  beginning  at  the  5-foot  position,  as  shown  in 
figure  2,  mean- velocity  distributions  were  determined  at  x  =  7*5,  8.5, 
9-5,  and  10.5  feet  at  a  free-stream  speed  of  108  feet  per  second.  These 
are  shown  in  figure  5-  A  lower  value  of  U-^  would  have  shifted  the 

transition  point  too  close  to  the  trailing  edge  of  the  plate.  Figure  6 
shows  the  same  velocity  profiles  on  a  nondimens ional  basis.  In  order 
to  avoid  the  arbitrariness  of  5,  the  nondimensional  distance  from  the 
wall  y/6*  is  vised.  It  is  seen  from  figure  6  that  downstream  from  the 
extreme  position  of  the  transition  point,  that  is,  x  =  6.75  feet",  the 
profiles  exhibit  similarity.  The  two-dimensional  nature  of  the  flow 
was  checked  by  making  velocity  measurements  10  inches  below  and  above 
the  center  line  at  x  =  10.5  feet  and  comparing  with  the  velocity  dis¬ 
tribution  at  the  center  line.  Agreement  among  the  three  sets  of  distri¬ 
butions  confirmed  the  two-dimensionality  of  the  flow. 

The  velocity  measurements  for  the  0.04-  and  0.25- inch  rods  are 
shown  in  figures  7  and  8.  At  a  free-stream  speed  Uq_  =  108  feet  per  sec¬ 
ond,  the  flow,  after  passing  over  the  rods,  became  turbulent  in  the 
separated  layer  and  reattached  to  the  surface  in  the  turbulent  state. 

At  a  lower  free-stream  speed  U-L  =  30  feet  per  second  for  the  0.04- inch 

rod,  visual  observation  of  an  oscilloscope  trace  obtained  by  traversing 
downstream  from  the  rod  at  a  distance  of  0.03  inch  from  the  surface,  with 
a  hot-wire  anemometer,  showed  the  reappearance  of  laminar-boundary- layer 
oscillations,  indicating  that  the  separated  layer  returned  to  the  surface 


v 


as  a  laminar  layer  without  the  occurrence  of  transition.  This  effect 
was  also  noticed  by  Liepmann  and  Fiia  (reference  10).  Figures  9  and  10 
give  some  of  the  representative  profiles  for  the  0.04-  and  0.25-inch 
rod  in  nondimensional  form.  The  profile  for  x  =  5 • 17  feet  in  figure  9 
differs  from  those  at  6.0  and  10.5  feet  by  more  than  the  experimental 
scatter.  This  would  place  the  beginning  of  similarity  for  the  0.04-inch 
rod  between  x  =  5.17  feet  and  x  =,6.0  feet.  In  figure  10  there  is  no 
evidence  of  the  beginning  of  similarity  downstream  of  the  0.25- inch  rod, 
and  comparison  of  the  profile  at  x  =  10.5  feet  with  other  distributions, 
at  the  same  Reynolds  nuniber,  which  do  exhibit  similarity,  indicated 
that,  even  at  x  =  10.5  feet,  6.5  feet  downstream  from  the  rod,  the 
velocity  distribution  is  not  yet  similar.  Figures  ll(a)  and  11(b)  show 
the  distributions  for  the  0.04-,  0.11-,  0 . l8- ,  and  0.25- inch  rods  at 
x  =  6.0  feet  and  x  =  10.5  feet,  respectively.  It  is  seen  that  the 
deviation  among  the  profiles  is  less  at  x  =  10.5  feet  than  at  6.0  feet, 
illustrating  the  trend  toward  similarity  and  the  dependence  of  the  simi¬ 
larity  position  on  the  rod  diameter.  The  trend  in  figure  10  is  in  the 
opposite  direction  of  a  Reynolds  number  effect,  but  the  Reynolds  number 
change  involved  is  small  enough  to  be  negligible  compared  with  the  dif¬ 
ferences  indicated.  It  is  expected  that  immediately  downstream  from 
the  rod  there  would  be  an  influence  on  dP/dx  and  dP/dy,  but  there 
is  nO  evidence  that  it  persisted  to  any  great  length  downstream.  If  the 
similarity  position  for  the  0.04- inch  rod  is  taken  at  x  =  5*5  feet, 
that  is,  1.5  feet  downstream  from  the  rod,  the  ratio  of  similarity  dis¬ 
tance  to  rod  diameter  is  450.  Assuming  a  linear  effect  of  the  rod 
diameter  on  the  similarity  distance  would  put  the  similarity  position 
for  the  0 . 11- ,  0 . l8- ,  and  0.25-inch  rods  at  x  =  8.2,  10. 9,  and 
13.5  feet,  respectively.  This  is  applied  only  for  the  condition  of  the 
rod  at  4  feet  and  U1  =  108  feet  per  second.  The  wind  speed  and  the 

position  of  the  rod  in  the  layer,  that  is,  the  size  of  the  rod  compared 
with  the  boundary- layer  thickness,  may  be  other  variables  which  have  an 
effect.  The  conclusion  from  this  is  that  too  much  length  is  required 
for  distortions  to  disappear  for  rods  to  be  practical  thickening  devices. 

The  velocity  distributions  at  a  free- stream  speed  =  43  feet 

per  second,  for  a  24-mesh  screen,  making  an  angle  of  0.5°  with  the 
plate,  are  shown  in  figure  12.  Although  no  summary  plot  is  given,  it 
may  be  stated  that  the  profiles  from  x  =  3*5  feet  to  x  =  10.5  feet 
exhibited  similarity.  The  boundary- layer  thickness  at  x  =  10.5  feet 
is  approximately  3  inches  as  compared  with  1.4  inches  for  the  condition 
of  free  transition.  The  effect  of  wider  angle  or  greater  solidity  was 
to  move  the  similarity  position  farther  downstream  and  to  give  a  thicker 
layer.  For  example,  the  profiles  produced  by  a  24-mesh  screen,  at  an 
angle  of  2.5°  with  the  plate,  did  not  show  similarity  at  x  =5*0  feet. 

At  x  =  9*5  feet  the  flow  was  similar,  and  the  layer  was  5*2  inches 
thick.  Since  no  measurements  were  made  between  x  =  5«0  feet  and 
x  =  9.5  feet,  the  position  at  which  similarity  first  occurred  is  unknown. 
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At  an  angle  of  5°  the  flow  was  no  longer  similar  for  the  length  of  the 
plate.  The  use  of  screens  seemed  to  he  an  effective  method  of  thicken¬ 
ing  hut  was  not  further  investigated  because  of  the  difficulties  entailed 
in  mounting  a  screen,  keeping  it  under  tension,  and  maintaining  its  shape 
for  long  periods  of  time. 

Figures  13 ;  l4,  and  15  show  the  velocity  distributions  downstream 
from  the  distributed  roughness  at  free-stream  speeds  of  35*  55;  and 
108  feet  per  second.  Since  the  general  characteristics  for  the  three 
speeds  are  similar,  the  nondimensional  profiles  are  given  in  figure  1 6 
only  for  =  55  feet  per  second.  It  is  seen  that  a  distance  of  about 
3-0  feet  was  required  for  similarity  and  that  at  x  =  10.5  feet  the 
layer  was  approximately  3  inches  thick,  which  for  =  55  feet  per  sec¬ 
ond  is  the  thickness  of  a  layer  having  its  beginning  about  ik  feet 
upstream.  In  other  words,  the  virtual  origin  of  the  layer  was  dis¬ 
placed  3*5  feet  upstream  of  the  leading  edge.  There  was  some  gain  in 
effective  distance  at  the  higher  speed.  At  =  108  feet  per  second, 

the  virtual  origin  was  about  5*5  feet  upstream  of  the  leading  edge. 

Comparison  of  the  velocity  distributions  in  the  various  artificially 
thickened  layers  is  made  with  that  of  free  transition,  in  figure  17;  at 
approximately  the  same  Reynolds  number  Rg.  The  agreement  is  excellent. 

To  test  the  universal  character  of  the  velocity  distributions,  it  is  of 
interest  to  compare  them  with  those  of  other  investigations.  In  fig¬ 
ure  18,  the  distributions  in  the  layer  thickened  by  sand  roughness  are 
compared  with  the  measurements  of  Hama  (reference  ll)  and  Schultz-Grunow 
(reference  12).  The  two  curves  show  the  effect  of  Reynolds  number,  and 
the  data  are  sufficiently  consistent  to  conclude  that  the  velocity  dis¬ 
tributions  obtained  for  the  artificially  thickened  layer  are  character¬ 
istic  of  a  turbulent  boundary  layer  with  zero  pressure  gradient  and  a 
smooth  surface. 


Comparison  with  Logarithmic  Laws 

Much  of  the  theoretical  work  that  has  been  done  on  the  turbulent 
boundary  layer  is  based  on  the  application  of  the  logarithmic  laws  of 
velocity  distribution  derived  for  flow  in  pipes  and  channels  (refer¬ 
ence  13).  These  involve  the  adoption  of  the  mixing-length  concept,  a 
form  for  the  shearing  stress  distribution,  and  the  omission  of  the 
influence  of  viscosity  on  the  turbulence. 

The  wall-proximity  law  and  the  velocity-defect  law  for  the  Velocity 
distribution  are  compared  with  experimental  results  in  figures  19  and  20. 
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The  value  of  UT  in  the  present  measurements  was  determined  from  the 

measured  velocity  profile  by  means  of  the  Von  Karmlin  momentum  equation 
which  for  zero  pressure  gradient  gives 


The  mean- velocity  distributions  for  the  artificially  thickened  turbu¬ 
lent  boundary  layer  with  zero  pressure  gradient,  together  with  those 
of  Schultz-Grunow  (reference  12)  and  Freeman  ( ref erence  l4) ,  are  plot¬ 
ted  in  the  form  of  U/UT  against  yUT/v  in  figure  19*  In  contrast 

with  the  measurements  for  pipes  and  channels,  it  is  seen  that  there  is 
a  systematic  increase  in  the  deviation  from  the  logarithmic  law  as  the 
Reynolds  number  decreases.  In  the  region  where  the  law  is  valid  fig¬ 
ure  19  gives  values  of  4.2  for  A  and  5*75  for  B.  The  following  table 
gives  the  comparison  with  the  measurements  in  smooth  pipes  and  channels 
for  wall  proximity. 


Investigator 

Nikuradse  (reference  15) 
Nikuradse  ( reference  1 6) 
Donch  (reference  17) 
Laufer  (reference  5) 
Klebanoff  and  Diehl 


Type 

A 

B 

Pipe 

5.0 

5-75 

Channel 

5.0 

!•  5-75 

Channel 

.  5.0 

5-75 

Channel 

7.2 

5-75 

Boundary- layer 

b.2 

5-75 

The  commonly  used  values  for  the  experiments  of  Nikuradse  in 
smooth  pipes  for  points  near  the  wall  is  5*5  for  B  and  5-8  for  A, 
but  Millikan  (reference  l8)  has  shown  that  values  of  5*75  and  5*0  are 
in  closer  agreement  with  the  measurements.  The  broken  line  in  figure  19 
shows  that  the  power  law 


=  8.1 


<?)’ 


is  a  good  approximation  for  the  higher  Reynolds  number  Rg  =  152,000 

throughout  most  of  the  layer,  but  it  should  be  remarked  that  at  the 
lower  Reynolds  numbers  the  agreement  was  not  so  good.  The  distribu¬ 
tions  at  lower  Reynolds  numbers,  on  a  log- log  basis,  had  a  sinuous 
appearance  which  seemed  to  decrease  with  increasing  Reynolds  numbers. 
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In  figure  20,  both  boundary- layer  and  channel  results  are  plotted 
in  the  form  suggested  by  the  velocity-defect  law.  It  is  seen  that 
V  results  for  both  the  boundary  layer  and  the  channel  fit  the  straight 

lines  up  to  y/8  =  0.15*  The  value  of  B  seems  to  be  universal  for 
boundary  layers  and  channels  within  this  region.  While  results  for 
*  pipes  are  omitted,  it  is  remarked  that  the  same  holds  true  for  pipe  flow. 

It  will  be  noted  that  the  observed  results  are  independent  of  Reynolds 
number  when  plotted  in  this  way,  indicating  separate  universal  relation¬ 
ships  of  the  form 

U-l  -  U 
UT 

A  point  to  be  made  with  emphasis  is  that  boundary- layer  profiles  are 
not  like  channel  and  pipe  flow  profiles,  and  it  is  seen  that  an  arti¬ 
ficially  thickened  boundary  layer  conforms  to  the  pattern  established 
in  figures  19  and  20  after  the  similarity  stage  has  been  reached. 


Spectrum  and  Intensity  of  u-Fluctuation 

It  has  been  shown  in  the  previous  section  that,  at  a  sufficient 
distance  downstream,  the  mean-velocity  distributions  assume  a  universal 
character.  However,  the  lack  of  knowledge  concerning  the  relationship 
between  the  turbulence  and  the  mean  flow  prohibits  the  use  of  mean- 
velocity  distributions  as  direct  evidence  that  the  turbulence  is  fully 
developed.  Since  it  is  impractical  to  measure  all  the  characteristics 
of  the  turbulence,  and  in  view  of  the  interest  in  measurements  of  the 
spectrum  in  shear  flow,  the  one-dimensional  energy  spectrum  was  chosen 
as  the  basis  of  comparison  for  the  turbulence  in  the  artificially 
thickened  layers  with  that  resulting  from  free  transition. 

The  mean-square  value  of  the  u-fluctuations  was  measured  between 
frequencies  n  and  n  +  dn  at  values  of  n  from  10  to  10,000  cycles 
per  second.  Over  most  of  the  frequency  range  the  accuracy  of  the  spec¬ 
trum  measurements  is  within  ±10  percent.  The  accuracy  was  somewhat  less 
at  the  low-frequency  end  below  100  cycles  per  second  because  of  large- 
amplitude  fluctuations  and  at  the  high  end  above  6000  cycles  per  second 
because  of  noise  and  possible  wire-length,  effect.  No  wire-length  cor¬ 
rections  were  made,  and  the  noise  level  at  10,000  cycles  per  second  was 
approximately  50  percent  of  the  total  signal.  A  direct  over-all  meas¬ 
urement  of  the  intensity  u’/l^  was  made  during  each  spectrum  determi¬ 
nation  to  compare  with  that  calculated  from  the  spectrum,  and  to  serve 
as  a  check  on  the  accuracy.  In  general,  the  agreement  was  within 
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'  5  percent.  The  spectrum  curves  are  given  in  the  n< 

against  Here  F(n)  is  the  normalized 

hx  u 

defined  by 

J1  F(n)  dn  =  1 


and  Lx  is  the  integral  scale  defined  by 


-X 


Lx  = 


where  Rx  is  the  longitudinal  space  correlation.  The  value  of  Lx  was 

obtained  from  the  intercept  of  the  F(n)  curve  at  n  =  0  on  the 
assumption  that  a  space  correlation  was  related  to  a  time  correlation 
by  x  =  Ut  and  using  the  Fourier  transform  written  in  the  following 
form: 


L*  =_ I - - 

The  values  of  Lx  calculated  by  equation  (3)  for  the  turbulent  boundary 
layer  resulting  from  free  transition  and  the  boundary  layers  artificially 
thickened  with  a  0.04- inch  rod  and  by  roughness  as  described  earlier  are 
given  in  figures  21,  22,  23,  and  24.  It  is  noted  that,  except  near  the 
wall.  Lx  is  constant  across  the  layer  and  is  proportional  to  the 
boundary- layer  thickness,  having  a  value  of  1^/5  of  approximately  0.4. 
It  should  be  pointed  out  that  the  exactness  of  equation  (2)  is  open  to 
question  at  high  turbulence  levels,  since  it  requires  that  the  fluctua¬ 
tions  at  a  point  be  the  result  of  a  pure  translation  with  velocity  XJ 
of  a  fixed  turbulence  pattern. 
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For  the  condition  of  free  transition,  figure  21  shows  a  set  of 
curves  at  x  =  10.5  feet,  with  =  108  feet  per  second,  at  different 

distances  from  the  surface.  It  is  seen  that  there  is  a  systematic 
variation  across  the  "boundary  layer  up  to  y/&  =  0.57.  Since  Lx  was  - 

used  as  the  reference  length,  all  curves  are  brought  together  at  zero 

*  frequency.  The  trend  in  going  toward  the  surface  is  for  the  higher 
frequencies  (smaller  eddies)  to  have  proportionately  more  of  the  energy. 

It  may  be  that  near  the  surface  some  of  the  smaller  eddies  can 
derive  energy  from  the  steep  velocity  gradient  which  exists  and  that, 
farther  away  from  the  plate  where  the  velocity  gradient  is  comparatively 
small,  local  isotropy  begins  at  a  lower  wave  number.  It  is  hoped  that 
future  measurements  of  the  spectrum  of  shearing  stress  will  throw  some 
light  on  this  point. 

Comparison  is  made  in  figure  22,  at  x  =  10.5  feet  and  y/8  =  0.03 
and  y/8  =  0.57,  with  the  spectrums  obtained  in  the  boundary  layer 
resulting  from  free  transition  and  those  obtained  in  the  artificially 
thickened  layers  using  a  0.04-inch  rod  and  sand  roughness.  The  free - 
stream  speed  was  108  feet  per  second  for  the  rod  and  free  transition 
and  55  feet  per  second  for  the  roughness  so  that  the  boundary-layer 
Reynolds  number  is  approximately  the  same.  The  spectrums  are  similar 

*  to  one  another  and  the  same  family  of  spectrums  exists  in  the  artifi¬ 
cially  thickened  layers  as  for  that  resulting  from  free  transition. 

*  Spectrum  measurements  made  at  various  distances  downstream  from 
the  0. Ok- inch  rod,  at  y/8  values  of  0.57  and  0.4  and  U  =  10.8  feet 

per  second,  are  given  in  figure  23-  It  will  be  noted  that  the  curves 
for  x  =  6.0,  7 • 5 ,  and  10.5  feet  are  similar.  Although  a  single  curve 
is  drawn  through  the  measured  points,  there  is  evidence  of  a  Reynolds 
number  effect  in  the  downstream  direction  on  the  high-frequency  end  of 
the  spectrum.  The  energy  associated  with  the  higher  frequencies  is  so  ' 
small  'that  it  does  not  distort  the  normalized  curve  at  the  low-frequency 
end.  The  curve  for  x  =  4.17  feet  at  y/8  =  0.4  which  is  only  2  inches 
from  the  rod  is  definitely  dissimilar.  This  curve  is  practically  flat 
up  to  1500  cycles  per  second.  The  same  characteristic  flatness  was  also 
observed  at  y/8  =  0.03,  so  no  particular  significance  should  be  attached 
to  the  value  of  y/8  =  0.4,  as  compared  with  the  rest  of  the  curves  with 
y/8  =  0.57.  Measurement  of  the  spectrum  for  the  condition  of  free  tran¬ 
sition,  at  x  =  6.9  feet,  which  was  as  close  as  it  was  possible  to 
approach  the  transition,  because  of  the  back  and  forth  movement  of  tran¬ 
sition,  showed  no  evidence  of  the  characteristic  flatness  observed  close 
to  the  rod.  In  fact,  it  gave  evidence  that  similarity  of  the  spectrums 
began  at  the  downstream  extreme  of  the  transition  region.  In  the  case 

**  of  free  transition,  the  spectrum  was  also  measured  in  the  laminar  layer 

at  x  =  4.5  feet.  A  very  sharp  peak  was  found  at  280  cycles  per  second. 
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the  predominant  frequency  of  the  laminar-boundary- layer  oscillations, 
and  little  energy  was  found  at  other  frequencies.  There  was  no  evidence 
of  this  peak  in  the  turbulent  boundary  layer  at  x  -  6.9  feet.  It  was 
thought  at  first  that  the  flatness  of  the  spectrum  observed  near  the 
rod  was  associated  with  the  distorted  mean- velocity  profile  and  that 
here  was  evidence  of  a  connection  between  turbulence  distortion  and 
mean- velocity  distortion.  However,  figure  24,  which  pertains  to  thick¬ 
ening  by  roughness,  shows  no  effect  in  the  spectrum  at  the  3-foot  posi¬ 
tion  where  the  mean-velocity  profile  was  definitely  distorted.  It  wou^d 
seem  that  the  flatness  of  the  spectrum  measured  near  the  rod  is  an  effect 
peculiar  to  the  rod  and  may  be  associated  with  eddy  frequencies  from 
the  rod  superimposed  on  a  more  or  less  normal  spectrum. 

It  may  be  expected  that,  in  the  turbulent  boundary  layer,  at  a 
distance  far  enough  away  from  the  wall,  an  extensive  region  of  local* 
isotropy  exists.  It  is  of  interest  to  see  how  well  the  spec trums  in 
the  artificially  thickened  turbulent  boundary  layer  wnere  the  Reynolds 
number  of  turbulence  is  high  agree  with  existing  theories  predicting 
the  shape  of  the  spectrum  in  isotropic  turbulence. 

Before  making  the  comparison  it  may  be  well  to  discuss  briefly  _ 
some  important  features  of  the  theory.  A  fundamental  concept  in  the 
theories  introduced  by  Kolmogoroff  (references  1  and  2)  and  Heisenberg 
(reference  19)  is  that  energy  passes  into  the  spectrum  by  way  of  the 
larger  eddies  and  is  then  passed  on  to  smaller  and  smaller  eddies. 
Appreciable  dissipation  by  viscosity  starts  only  when  the  eddies  get 
small  enough.  Beyond  a  certain  wave  number,  smaller  eddies  get  their 
energy  only  from  larger  eddies  and  are  supplied  sufficient  energy  to  be 
kept  in  statistical  equilibrium  even  in  cases  where  the  larger  eddies  are 
decaying.  When  the  Reynolds  number  of  turbulence  u*L xjv  is  high,  there 
will  exist  in  the  lower  wave-number  part  of  the  equilibrium  range  a  sub¬ 
range  in  which  inertial  forces  predominate  and  the  transfer  of  energy 
is  equal  to  the  total  dissipation  e,  and  at  the  high  wave-number  end, 
another  subrange  where  viscous  effects  predominate.  Heisenberg  assumed 
that  the  transfer  of  energy  at  wave  number  k  was  caused  by  a  turbu¬ 
lence  friction  produced  by  eddies  with  wave  numbers  greater  than  k  and 
that  the  total  dissipation  involved  a  turbulence  viscosity  plus  ordinary 
viscosity  in  an  expression  of  the  following  form: 
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Here  the  second  term  within  the  brackets  represents  the  turbulence  vis- 
,  cosity,  and  x  is  an  absolute  constant.  The  solution  of  equation  (4) 
gives  the  spectrum  in  the  equilibrium  range  the  form 


where  ks  is  a  wave  number  in  the  intermediate  range  where  inertial 
forces  and  viscosity  both  play  a  role,  and  is  given  by 


.  The  equilibrium  spectrum  given  by  equation  (5)  deals  with  a  three- 

dimensional  energy  spectrum.  In  order  to  compare  with  experiment  it  is 
necessary  to  transform  this  equation  to  the  one-dimensional  form  using 
i  the  following  relationship  also  derived  by  Heisenberg: 


The  one-dimensional  form  of  equation  (5)  is  found  to  be 
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where 


arid 

ksx  =  0.645  ks 

In  the  inertial  range,  that  is,  kx  «  ksx,  equation  (7)  becomes 


Ex(kx)  =  jk  A2k* 


(9) 


where 


In  figure  24,  equations  (8)  and  (9)  are  compared  with  the  spectrums 
measured  at  different  distances  downstream  from  the  roughness,  at 
y/S  =0.57  and  1^=  55  feet  per  second.  At  the  high- turbulence 
Reynolds  number  of  the  boundary  layer,  the  existence  of  a  -5/3  region 

is  quite  evident.  The  agreement  with  the  n"^  law  should  be  taken  with 
some  reservation  because  of  the  lesser  accuracy  of  the  measurements  at 
the  high-frequency  end  of  the  spectrum.  It  should  be  mentioned  in  pass¬ 
ing  that  the  slope  of  the  spectrum  curves  at  y/&  =  0.03,  such  as  shown 
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4 


in  figure  22,  is  approximately  one,  where  the  -5/3  region  exists  at 
y/8  =  0.57. 


It  is  seen  from  figure  25  that  in  the  intermediate  range  of  fre¬ 
quencies  equation  (7)  does  not  agree  so  well  as  the  interpolation 
formula  ' 


suggested  by  Heisenberg.  For  convenience,  equations  (7)  and  (10)  are 
represented  in  figure  25  as 


and 


(10a) 


where 


F(z)  - 

Lx 


and 


nLx 
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and  A0  and  zs  are  experimentally  determined.  Formula  (10a),  as 

shown  in  figure  2 6,  is  also  in  good  agreement  with  some  unpublished 
measurements  made  in  isotropic  turbulence  50  inches  downstream  from  a  ^ 

1-inch  grid  in  the  KBS  -foot  wind  tunnel.  The  fact  that  the  analyt- 

d  * 

ical  relation  (7)  does  not  agree  with  the  measured  spectrums  so  well  as 
the  empirical  expression  (10)  would  indicate  that  equation  (k)  may  not 
yet  be  quite  exact.  Expressions  (7a)  and  (10a)  are  identical  in  the 
extreme  regions  of  the  equilibrium  range,  and  it  is  of  interest  to  cal¬ 
culate  the  value  of  x  from  the  constants  A0  and  zs.  Since 

Aq  and  zg  are  determined  from  a  nondimens ional  normalized  represen¬ 
tation  they  are  connected  to  X  by  the  following  relations. 

The  term  E(k)  of  the  three-dimensional  spectrum  is  defined  such 

that 
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since 


kv  = 


2jm  2wz 
U  Lx 


4  2n 


gF(z)  =Ex(kx) 


In  the  inertial  range. 


F(z) 


AqzS 


5/3z-5/3 


^x(^x)  ^ozs 


5/3 


i 


and 


110 

9 


AqZs 


5/3 


(11) 


It  will  be  noted  that  e  and  '  u2  must  be  determined  in  addition 
to  the  quantities  already  mentioned  before  X  can  be  calculated  by 
equation  (ll).  For  isotropic  turbulence,  the  expression  for  the  dis¬ 
sipation  (reference  20)  i °- 

€  =  15v4  (12) 

X d 

where  X  is  the  microscale,  obtainable  from  a  spectrum  measurement  by 


a 


n2F(n)  dn 


(13) 


Partial  justification  for  the  use  of  such  relationships  where  the 
turbulence  does  not  have  general  isotropy  is  afforded  by  figure  27  > 
where  it  is  seen  that  about  90  percent  of  the  integral  in  equation  ( 13) 
is  contributed  by  frequencies  above  500  cycles  per  second.  This  is  in 
bhe  region  where  the  -5/3-  and  -7-power  laws  hold  and  is  therefore  the 
region  of  local  isotropy.  It  is  concluded  therefore  as  a  reasonable 
approximation  that  X  determined  by  equation  { 13)  is  the  appropriate 
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and  only  scale  entering  into  the  viscous  dissipation.  However,  much  of 
the  value  of  u2  is  contributed  below  500  cycles  per  second,  and  the 

assumption  that  u2  +  v2  +  w2  equals  3u2,  as  made  in  equation  (12), 

is  not  a  sufficiently  good  approximation.  Since  v2  and  w^  were  not 
measured,  the  results  of  reference  (6)  were  used  to  obtain 

u2  +  v2  +  w2  =  0.62  (3u2) 

at  y/B  =  0.57.  The  factor  0.62  is  applied  to  the  right  side  of  equa¬ 
tion  (12)  and  to  equation  (ll)  where  u2  appears. 

The  value  of  X  calculated  from  equation  (ll)  for  the  turbulent 
boundary  layer  artificially  thickened  by  sand  roughness,  at  x  =  10.5  feet 
and  y/8  =  0.57,  with  the  values  of  A0  and  zs  given  in  figure  26, 
is  0.23.  The  value  of  X  calculated  for  the  case  of  isotropic  turbu¬ 
lence  shown  in  figure  26  is  O.65.  By  a  similar  procedure  X  was  cal¬ 
culated  for  the  spectrum  measured  by  Laufer  in  fully  developed  channel 
flow  (reference  5) ,  and  a  value  of  0.52  was  obtained.  At  present  there 
are  insufficient  experimental  data  to  determine  the  reason  for  the  var¬ 
iations  obtained.  It  is  of  interest  to  mention  that  the  turbulence 
u*l 

Reynolds  numbers  —  -x ■  were  1790 ,  550,  and  190  for  the  boundary  layer, 

channel,  and  isotropic  turbulence,  respectively,  and  that  corresponding 
increasing  values  of  ,X  were  obtained; 

The  insensitivity  of  the  spectrum  to  distortions  in  the  mean- 
velocity  distributions  as  illustrated  by  figure  2k  made  it  desirable 
to  check  further  on  the  fully  developed  nature  of  the  turbulence  in  the 
artificially  thickened  layer  by  measuring  the  distribution  of  u'/U]_. 

Intensity  measurements  were  made  at  various  distances  downstream  from 
the  roughness  at  =  108  feet  per  second,  and  these  are  given  in  fig¬ 
ure  28.  At  x  =  3.0  feet,  1  foot  downstream  from  the  edge  of  the  rough¬ 
ness,  the  intensity  is  flat  up  to  y/B  =  0.4  and  is  apparently  the 
result  of  the  diffusion  of  the  high  turbulence  generated  by  the  rough¬ 
ness.  From  x  =  5.0  to  x  =  10.5  feet,  the  distributions  of  uJ 
exhibit  similarity,  which  is  the  same  region  for  which  the  mean- velocity 
distributions  were  similar.  The  change  in  turbulence  level  accompanying 
a  change  in  speed  from  108  to  55  feet  per  second  is  also  shown.  It 
should  be  noted  that  there  is  no  change  in  level  accompanying  the  vari¬ 
ations  in  boundary- layer  Reynolds  number  due  to  changes  in  boundary- 
layer  thickness. 
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An  interesting  feature  of  the  turbulence  in  the  boundary'  layer  was 
its  intermittent  nature  in  the  outer  portions  of  the  layer.  This  was 
also  observed  by  Townsend  for  a  wake  (reference  .21)  and  by  Corrs'in  for 
a  turbulent  jet  (reference  22).  It  seems  to  be  characteristic  of  shear 
flow  with  a  free  boundary.  Figure  29  shows  two  oscillograms  taken  at 
x  =  10.5  feet,  at  y/S  =0.03  and  y/5  =0.75-  At  y/S  =0.03  there 
is  no  trace  of  intermittency,  while  at  y/8  =  0.75  it  is  quite  marked. 
Observation  of  the  oscilloscope  trace  indicated  there  is  a  gradual 
increase  in  the  intermittency  as  the  distance  from  the  surface  increased. 
It  is  planned  in  the  future  to  make  some  quantitative  measurements  of 
the  intermittency,  and  its  effect  on  the  measured  turbulence  intensity. 


CONCLUSIONS 


From  tests  made  to  determine  the  feasibility  of  artificially 
thickening  a  turbulent  boundary  layer  on  a  flat  plate,  the  following 
conclusions  can  be  drawn:  . ,  . 

1.  It  has  been  established  that  by  artificial  methods  it  is  pos¬ 
sible  to  thicken  substantially  a  fully  developed  turbulent  boundary 
layer  on  a  flat  plate.  The  layer  so  formed  becomes  characteristic  of 
that  developed  from  the  beginning  on  a  smooth  flat  plate  after  passing 
over  the  required  length  of  smooth  surface.  For  a  given  boundary- layer 
thickness,  artificial  thickening  makes  possible  a  saving  in  length 
provided  that  there  is  enough  length  for  the  layer  to  become  free  from 
distortions. 

2.  Almost  any  kind  of  obstruction  placed  on  the  surface  would 
produce  a  thickened  turbulent  flow,  but  only  from  a  certain  class  did 
downstream  effects  disappear  in  a  reasonably  short  distance.  This  class 
appeared  to  be  characterized  by  the  ability  to  absorb  energy  from  the 
mean  flow  without  introducing  large-scale  disturbance.  Two  examples  of 
this  class  are  as  follows: 

(1)  A  fine-mesh  screen  placed  so  that  the  plate  and  screen 

form  a  narrow  wedge  opening  into  the  wind 

(2)  Sand  roughness  cemented  to  the  surface 

3.  The  mean- velocity  distribution  appears  to  afford  a  simple  means 
of  identifying  a  fully  developed  turbulent  boundary  layer. 

4.  The  logarithmic  laws  of  velocity  distribution  are  not  in  good 
agreement  with  the  measured  distributions  in  the  turbulent  boundary 
layer  with  zero  pressure  gradient  and  are  applicable  only  in  a  small 
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region  near  the  wall.  The  velocity  profiles  in  the  boundary  layer  are 
independent  of  Reynolds  number  on  a  velocity-defect  basis,  indicating 
a  universal  relationship  of  the  form  « 

as-fww 

UT 

where  U  is  the  mean  velocity  in  the  boundary  layer,  is  the  mean 

velocity  in  the  free  stream  just  outside  the  boundary  layer,  UT  is 

the  square  root  of  the  ratio  of  the  shearing  stress  at  the  surface  to 
the  density,  y  is  the  distance  normal  to  the  surface  measured  from 
the  surface  of  a  flat  plate,  and  6  is  the  boundary-layer  thickness. 

I 

5.  The  longitudinal  scale  of  turbulence  Lx  in  a  fully  developed 

layer  is  proportional  to  the  boundary- layer  thickness  5,  with  a  value 
of  Lx/8  of  approximately  0.4,  and  except  near  the  wall  is  constant 

across  the  layer. 

6.  There  is  a  systematic  variation  in  the  one-dimensional  energy 
spectrum  across  the  boundary  layer  up  to  y/S  =  0.6,  and  in  the  outer 
region  of  the  boundary  layer,  that  is,  y/&  >  0.6,  the  existing  theories 
for  the  spectrum  in  isotropic  turbulence  can  be  applied. 

*  7.  There  is  an  intermittent  character  to  the  flow  in  the  turbulent 

boundary  layer  similar  to  that  found  for  wakes  and  jets. 

National  Bureau  of  Standards 

Washington,  D.  C.,  September  15,  1950 
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(a)  Pressure  distribution  along  flat  plate. 


(b)  Variation  of  dynamic  pressure  as  measured  with  a  surface  tube  through 
transition  region,  =  108  feet  per  second. 

Figure  2.-  Pressure— distribution  and  transition  data. 
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Figure  7.-  Mean-velocity  distributions  for  O.OU-inch  cylindrical  rod  in 
contact  with  surface  h  feet  from  leading  edge.  U]_  =  108  feet  per 


second. 


Figure  9.~  Nondimensional  mean-velocity  distributions  for  O.Ol^-inch 
cylindrical  rod  in  contact  with  surface  ij.  feet  from  leading  edge. 
U]_  =  108  feet  per  second. 
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Figure  10.-  Nondimensional  mean-velocity  distributions  for  0.25~inch 
cylindrical  rod  in  contact  with  surface  U  feet  from  leading  edge. 
Ui  =  108  feet  per  second. 
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Figure  12.-  Mean-velocity  distributions  for  1-foot  span  of  2l).-mesh 
screen  making  angle  of  0.5°  with  plate.  U]_  =  U3  feet  per  second. 
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Figure  1$.-  Mean-velocity  distributions  for  No.  16  sandpaper  covering 
first  2  feet  of  surface.  Up  =  108  feet  per  second. 


* 


Figure  17.-  Comparison  of  nondimensional  mean-velocity  distributions  having 

approximately  same  values  of  Rq. 
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Figure  19.-  Comparison  of  experimental  results  with  l/7-power  law  and 

logarithmic  law. 


velocity-defect  law, 


Figure  21. -  Spectrum  of  u-fluctuations  in  turbulent  boundary  layer 
resulting  from  free  transition,  x  =  10.5  feet;  Up  =  108  feet  per 
second.  . 
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Figure  22.—  Spectrum  of  u— fluctuations  in  turbulent  boundary  layer  at 
x  =  10.5  feet.  Ui  =  108  feet  per  second  for  free  transition  and 
O.OU-inch  rod;  Ui  =  55  feet  per  second  for  roughness. 
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Figure;  23.-  Spectrum  of  u-fluctuations  in  turbulent  boundary  layer  for’ 
O.OU-inch  rod  in  contact  with  surface  2;  feet  from  leading  edge. 

Up  =  108  feet  per  second.  V. 
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Figure  2k. -  Test  of  -5/3  law  and  -7  law  for  spectrum  of  u-fluctuation, 
in  turbulent  boundary  layer  thickened  by  roughness  for  first  2  fee 
of  surface.  U]_  =  55  feet  per  second. 


Figure  2 Comparison  of  measured  spectrum  of  u-fluctuations  with 
Heisenberg's  theoretical  formula  (broken  line)  and  interpolation  . 
formula  (solid  line).  Boundary  layer  thickened  by  No.  16  sandpaper 
for  first  2  feet  of  surface,  x  °  10.5  feet;  Ui  =  55  feet  per  second. 
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Figure  26.—  Spectrum  of  isotropic  turbulence  compared  with  Heisenberg's 
interpolation  formula.  Measurements  made  £0  inches  from  1-inch  grid 
in  uniform  stream. 
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Figure  29#-  Oscillograms  of  u-fluctuations  in  turbulent  boundary  layer, 
x  =  10.5  feet;  =  55  feet  per  second;  6  =  3.1  inches.  Time  interval 

between  dots,  l/60  second.  Boundary  layer  thickened  by  No.  16  sandpaper 
over  first  2  feet  of  surface. 
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